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We solve the Edwards- Anderson model (EA) in different Husimi lattices using the cavity method 
at replica symmetric (RS) and 1-step of replica symmetry breaking (1RSB) levels. We show that, 
at T = 0, the structure of the solution space depends on the parity of the loop sizes. Husimi lattices 
with odd loop sizes may have a trivial paramagnetic solution thermodynamically relevant for highly 
frustrated systems while, in Husimi lattices with even loop sizes, this solution is absent. The range 
of stability under 1RSB perturbations of this and other RS solutions is computed analytically (when 
possible) or numerically. We also study the transition from 1RSB solutions to paramagnetic and 
ferromagnetic RS solutions. Finally we compare the solutions of the EA model in Husimi lattices 
with that on the (short loops free) Bethe lattices, showing that already for loop sizes of order 8 both 
models behave similarly. 
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I. INTRODUCTION 



Spin glasses are among the most complex problems in Statistical Mechanics. During the 80's a lot of effort was 
devoted to the subject, see for example [1| for a comprehensive collection of relevant works during this decade. It was 
soon realized that the difficulties in finding an analytical solution to these problems depend strongly on the topology 
of the interactions between the variables [l], For example, in fully connected systems, in which each variable 
interacts with all the others, a compact solution may be found using the Parisi ansatzjl], 0,0). On the other hand, the 
situation for finite connectivity (FC) systems is far more complicated. The main difficulty is the appearance, after a 
standard replica calculation of an infinite number of over laps Hi]. Therefore, for many years only Replica Symmetric 
or variational solutions were known for these models. 

The importance of finite connectivity systems is twofold: first, one may hope to get a better understanding of 
finite dimensional systems, since (FC) models include the notion of neighborhood, a concept that is absent in fully- 
00 ] connected systems. Second, there is a clear connection between finite connectivity systems and many constraint 
satisfaction problems. For example, the K-sat @, the coloring Q, the traveling salesman 0] and the vertex cover 
' [Io| problems turn out to have a finite connectivity structure. 
qq . A few years ago, Mezard and Parisi[ll], [l2j generalized a technique already known as the cavity method[l3[ to 
deal with systems with many pure states. This generalization permitted, for the first time in FC systems, the formal 
introduction of replica symmetry breaking at different levels. Although, it is worth reminding that even the one-Step 
■ Replica Symmetry Breaking solution (1RSB) involves as an order parameter a functional distribution. 

Moreover, thanks to the aforementioned strong connection between finite connectivity spin systems and many 
constraint-satisfaction problems |14[ this approach sheds some light on the characteristics of the solution space of 
some of these problems [f| 0, M. llCf. (see also [HI, for recent developments in this field). In addition, this 
cavity method, inspired a novel message-passing algorithm to deal with single instances of several combinatorial 
problems [T3, El • Unfortunately, the occurrence of short loops in finite connectivity graphs introduces strong 
correlations among neighboring sites and in this case, the hypothesis behind the cavity method may be violated. As 
a consequence, message-passing algorithms usually fail when short loops are present [20l I2H I22I [23j . On the other 
hand, fully understanding the role of short loops in the energy landscape of finite dimensional spin glasses remains 
an elusive task[2^.[25j. 

To gain some insight about these problems, we think that it is convenient to look at the properties of a spin glass 
model, where the hypothesis behind the cavity approximation remain valid, but where the influence of short loop 
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structures may be analyzed in detail. Then, we choose to study the ground state characteristics of the Edwards- 
Anderson model in a Husimi graph. The tree-like structure of the Husimi lattice allows the use of the machinery 
behind the cavity method, and at the same time, we may tune the loop sizes and the connectivity of the graph to 
discover their influence on the zero temperature energy landscape of the model. 

The paper is organized as follows: In section [II] we introduce the model, the self-consistent cavity equations that 
solve it and the definition of the relevant physical quantities. Then, in section IIIII we show closed analytical results 
for the triangular Husimi lattice. In section HVl appears the analysis of more general lattices and finally in section fVl 
we present the conclusions of the work. 



A Husimi tree is formally "a connected graph where no bond (edge) lies in more than one cycle" . It can be visualized 
as a tree made out of loops, as shown in figure [TJ A Husimi tree is called pure if all the loops have the same length. 
It is also called regular if all the vertexes belong to the same number of loops [26| . On the other hand, the Husimi 
lattice (also Husimi graph) is a graph that looks locally like a Husimi tree, but where large loops are present. It is a 
random hyper-graph in the sense that by taking the short loops of the Husimi tree as building units, we define the 
ensemble of pure and regular Husimi graphs as the ensemble of all graphs in which the vertexes belong to K + 1 short 
loops of the same length c + 1 . 

For example, the simplest (and trivial) Husimi tree is the one in which c = 1 and coincides with a Cayley tree. 
The simplest Husimi lattice, c = 1 coincides with the Bethe lattice defined in [11] . The simplest, non-trivial, case 
of a Husimi lattice has c — 2, (a triangle) and K — 1, each vertex is shared by two triangles (see figure [T]). The 
generalization to more complex structures is straightforward, for example, either c or K or both may be taken as 
random variables such that the local structure of the graph changes from site to site. Here we consider the ensemble 
of pure and regular Husimi graphs, meaning that all short loops have the same length, and that all vertexes belong 
to the same number of (short) loops. 



FIG. 1: Left: schematic representation of a representative part of a Husimi graph with c = 2(triangle) and K = 1. Right: 
factor graph representation of the same Husimi lattice. In the factor graph representation, c + 1 is the degree of factor nodes 
(white squares) , while K + 1 is the degree of nodes (black circles) . 

In this work we study, using the cavity method at zero temperature [l2j , the ground state properties of the Edwards- 
Anderson model in an average Husimi graph. The Edwards- Anderson model is defined by the Hamiltonian: 



where <i,j> stands for the nearest neighbors in the graph, Si = ±1 are Ising variables located at the vertexes, and 
the bonds represent random exchange couplings Jjj. The couplings are taken from the distribution: 



II. THE MODEL AND THE CAVITY SOLUTION 
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where p € [— 1, 1] parameterizes the bias between ferromagnetic and anti-ferromagnetic interactions. If p = 1(— 1) the 
system is purely ferromagnetic (anti- ferromagnetic) , and if p = the system is unbiased. 

As will be discussed in more detail below, the cavity method may be easily written in term of messages between 
Husimi loops and sites. The idea is to use a factor graph that considers all the interactions around a loop as a single 
function node (see figure [T]) , then inheriting the tree- like structure of the Husimi graph at the loop level. We will 
use the cavity method at two levels of approximations: the replica symmetric approximation (RS) and the one step 
replica symmetry breaking approximation (1RSB). It will turn out that depending on c, K and p, either the RS 
approximation or the 1RSB are more appropriate to describe the ground state of the system. We leave for future 
works the analysis of the stability of the 1RSB solution. 



A. The Self-Consistent Cavity Equations at T = 

In order to write the cavity equations for the model, it is convenient to define the problem using a factor graph 
representation. A factor graph [27| is a bipartite graph in which one subset of the nodes represents local interactions 
while the other subset represents variables. In this case, it is reasonable to assume that each loop, a, corresponds to 
a function node, and each site, i to a variable node (see figure [1}. With this notation, equation (JTJ) may be written as 
a sum of loop contributions TL = £ 7 2 a H a where: 

c+l 

H a = — y^Jj t i+iSjSi-\-i (3) 

i 

where we assumed periodic boundary conditions in the index i within the loop. 

Then, since the resulting bipartite graph shares the structure of the Husimi lattice, it is locally tree-like and one 
expects the convergence of the cavity equations below. 



1. Replica Symmetric equations 

The self-consistent cavity equations in the Husimi lattice may be derived following the approach in [jjj]. Let us 
consider a set of c cavity spins S\ . . . S c each of which belongs to K loops and receives a field hi from the interactions 
with these loops. The iteration procedure consists in adding a new spin So to the cavity graph together with c+l 
couplings Jt,i+i, forming a loop (a function node) between all the c cavity spins and the spin So- With this new loop, 
the cavity spins S\ . . . S c complete their K + 1 loops, and a field contribution u a ^0i that can be computed from the 
values of h\ ... h c and Jq.i ■ ■ ■ J c .o, appears on spin So- 



u a ^ = u{J,h) = ^So min H a 
So Sl '" Sc 

1 c_1 
= ~- ^ S max { Jq.iSqSi + J lii+ iS l S l+ i + hjSj) + h c S c + J c .qS c S } (4) 
Z S\...S C * 

Then, the cavity field in the spin Sq is ho = Y]„—i k Ua wnere an u a are computed independently by the cavity 
iteration described above. The u are called cavity messages or biases, and can be interpreted as the contribution of 
a given loop a to the cavity field in So- The cavity fields hi, h2, ■ ■ ■ h c entering a loop, will be referred from now on 
as h for notation clarity. Similarly J will stand for the set of all the c+l coupling constants, J%j, of the function 
node. Equation Q is nothing but the difference between the minimum energy configuration of the spins at the loop 
when the cavity spin So is up and down. Following equation ^ it is easy to prove that independently of c and K the 
cavity messages in this model take only integer values: u 6 {—2, —1, 0, 1, 2}. 

Furthermore, it is expected that the distribution of these fields is stable during the iteration procedure [12j. This 
consideration yields the self-consistent cavity equations for the fields distribution V(h) and the messages distribution 
Q(u): 

K 

vih) = I s(h-j2 u «) n (5) 

a a—1 

Q{u) = Ej J S(u-u(J,h))f[dV(hi) (6) 
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V(h) and Q(u) are distributions of fields and biases that represent the probability of finding a field or bias in any site 
of the graph and Ej is an average over the Ji j. These two equations can be folded into one that is equivalent to the 
"Parent to Child" messages-passing algorithm described in [20] for single instances: 

K cK 

Q(u) = Ej / S(u - u{J, J2 «)) II dQ W ( ? ) 

"* »=i 

This is a self-consistent equation for the RS order parameter Q(u). As the biases are integers u G {—2, — 1, 0, 1, 2}, 
Q(u) can be parameterized with five numbers p-2,P-iiPo,Pi,P2- 



Q(u) = p- 2 5 u -2 + P-iSu-i + PoS U fl + pxS Ut i + p 2 S u ,2 (8) 

(or four if we consider the normalization constraint) that represent the probabilities p u of finding a message u going 
from a function node to a site. 

With the help of (JSJ) the self-consistent equation may be written as a set of five equations relating the parameters 
P_ 2 ,P-i,P0jP1jP2- Then, the right-hand side of (J7J) constitutes a sum of terms of degree c • K in the probabilities p, 
and degree c + 1 in p. These equations are highly coupled and non-linear: all 5 cK combinations of fields and 2 C+1 
combinations of couplings have to be considered. While in some cases one is able to find analytical solutions for all 
p, in more general situations equation (0) must be solved numerically. 

Once the RS solution Qrs(u) of the self-consistent equation flT]) for a given p is known, this can be used to compute 
the expected value of the energy in the graph, as: 

c/=-|±i Ej J a(j,f>)ndQK)+^^ / iE^in dfi K) ( q ) 

where the function a(-) is an energetic term very similar to the function u(-): 



a(J_,h) = — - max { J .i5 5i + Y^(J il j + iS',Si + i + hiSi) + h c S c + J c .oS c Sq} (10) 
Z Si—S c *r~t 

and is exactly the energy gained by the local state of the graph when the spin Sq is added to the cavity in the iteration 
procedure. 

The RS solution is valid if the system has a single pure state. In the more general case in which the system has 
many pure states it usually fails. The reason is that within the RS approximation one assumes that under the process 
of iteration the ground states of the graphs with TV and N + 1 spins are related. This is not necessary true[12J, and 
in such cases one must go beyond the RS solution to the 1RSB approximation. Another, perhaps, more numerical 
intuition is that running the message-passing algorithm, i.e. equation l[7]). over a single graph one observes that the 
RS equations have multiple fixed points, which is a signature of the appearance of multiple stationary states for the 
system (l6j. 



2. 1RSB equations 

The 1RSB cavity method assumes that there exists an exponential number A/*(e) ~ exp(NY,(e)) of pure states (ll. [l2| 
with intensive energy e = E a /N, characterized by the complexity function E(e). In this case the correct description 
of the interaction between a loop (node) and a spin is not given by a single message u, but by a distribution Q(u) 
representing all different messages in all pure states, and the correct order parameter is now the functional Q[Q(u)] 
defining the probability of finding a distribution Q(u) in a randomly selected node of the graph. The hypotheses made 
for the RS approximation are valid in each pure state a, but the statistics are more subtle, since the ground state can 
move from one pure state to another when the iteration procedure is carried out [ill . [l2j • The self-consistent equation 
for the new order parameter is similar to the RS equation for Q(u): 

. cK 

Q{Q] = Ej / 6 {F) (Q - QiLQj})) J[ dQ[Q a ] (11) 

"* a=l 
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where Q a represents the set of all the cK distributions Q a in the cavity iteration procedure. The functional delta 
accounts for the iteration of messages, and is defined by: 

i /■ K cK 

Q [l,Qa\ (u) = Jy-Qj J e^^>S(u - uU,J2u)) n dQ a (u a ) (12) 

where the exponential Boltzmann factor is the so called "re-weighting" that takes into account the energy shifts 
between states when the iteration procedure is carried out (T^ |. Note that the function d(-) was already defined in eq. 

m- 

In short, in the 1RSB approximation each of the states appears with a weight proportional to e M ^ where this fx 
plays the role of an inverse temperature and / is the corresponding free energy of the state [ill]. In the replica language 
the parameter fi stands for the product of the Parisi's replica parameter m and the inverse temperature (3 in the zero 
temperature limit (TH [l2l] . The term A[J_, Q a ] is a normalization constant given by: 



. cK 

A\J,Qa]= e^^l[dQ a (u) (13) 

J a 

and can be thought as a local partition function. 

The order parameter Qirsb[Q] solution of equation pip can be rarely found analytically, and furthermore, it 
depends on the variational parameter fi. The 1RSB solution, then, is quite more involved than the RS one, not only 
because the order parameter becomes a functional, but also because this functional has to be found numerically, and 
extremized over the parameter /x, as is explained next. 



3. Free energy, energy and complexity 

Within the 1RSB cavity formalism [ll], [l2j], a generalized free energy, 0(/i), can be defined as : 

(K + l)c 

<£(/x) =< <t>site(lj) > ~ — r~ < <t>node{v) > (14) 

C + l 

where the 4>node{fJ-) is the contribution from the loops and <p S it e (n) is the contribution coming from individual sites in 
the Husimi graph. The previous expression for 4>(fJ.) reduces to the one in [l2j for the EA model in the Bethe lattice, 
when c = 1. The site and node contributions are given by: 

c+l 

Ma{J,h)+\ho+u(l,h)\) JJdp.^.) ( 15 ) 

where Z S i te and Z noc i e are partition functions that assure the normalization of the probability distributions at site 
and loop levels respectively. 

Local distributions of fields P(h) and messages Q[u) hold the same relation than V(h) and Q(u) in equation 
though they have different meanings. Using this relation and equation (|13p . the computation of (j>((J,) in (TTi|) can be 
written in terms of Q[Q(u)] as: 

-Mm) = < io g ^[j,Q a ] > J Kc - V < i og />ii>«i TT dQiim) > (16) 

c+l c+l J - LJ - 

where < • > stands for the average over the order parameter distribution J Q[Q] as well as over the quenched disorder 
distribution Ej. This is not, however, a variational expression for and gives only valid results when the order 
parameter Q[Q] solves the self-consistent equation (fTTj) . The results obtained for the Replica Symmetric case can be 
recovered from the 1RSB expressions by taking the /i — > limit. 

The complexity is related to the free energy <j)(ix) via a Legendre transform [Til [T^] : 

S( M ) = n [e(/x) ~ 0(/i)] (17) 
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This justifies calling <fi a free energy, and £ and /x can be thought as the entropy and the temperature of a system 
whose configuration space is given only by the pure states of our original problem. On the other hand, given the relation 
of the complexity with the exponential abundance of states Af(e) ~ exp(NT,(e)), it is clear that the ground state of the 
system is one with zero complexity £ = 0. This defines the point p,* such that £(/i*) = and Uirsb = e(M*) = <Km*) 
is the ground state prediction at the 1RSB approximation. Furthermore, given the usual Legendre relations 



e(M) = WQi)] £( M ) = fjL^in) (18) 

it is clear that the point (j,* extremizes 4>{p). Following the analogy with the replica method P, H3 |. it can be 
shown that fi* actually maximizes <f>(n). As usual, it is sufficient to know <p(fi) to obtain all other thermodynamic 
potentials, but let us mention that explicit expressions for e(/x) and £(/x) in terms of the order parameter can be 
derived straightforwardly [III ■ 

To finish the presentation of the 1RSB cavity method some words must be said about the numerical method used to 
solve the equation (fTTj) . Since Q[Q(u)] is a mathematical obj ect very hard to deal with, this self-consistent equation 
is rarely solved analytically. The usual approach [Til . fl2L Il6l \l§L [28j is to represent the order parameter by a large 
population of distributions Q(u), and to use a fixed-point method to solve eq. (|11[) . In this work, each distribution 
Q{u) is represented by the five numbers p_2 ■ ■ -p 2 (similar to |J5])), and the order parameter consists of a population 
of N ~ 10 4 of such distributions. For each value of p, the solution of the self-consistent equation is found by replacing 
several times (~ 10 2 iV) a randomly selected member of the distribution by the result of ([T2"jl . After convergence, the 
free energy (flB)) is computed using other ~ 5 * 10 3 7V steps of the population dynamics algorithm. 



B. RS-1RSB Stability Analysis 

The stability under 1RSB perturbations of the RS solutions as a function of p can be studied with the method 
applied in [28l.l29|. The replica symmetric self-consistent equation can be recovered from the 1RSB [28], by restricting 
the space of the 1RSB order parameter: 

2 

Q[Q] - Qrs[Q] = £ P g S F (Q(u) - S u , q ) (19) 

q=-2 

In such a case, each distribution is forced to be deltaic and the re-weighting terms in (fT2|) can be factorized out of 
the integrals. Then, all 1RSB expressions turn to be exactly the ones obtained at the RS level. 

It is now fruitful to study the effect of the iteration procedure over a 1RSB parameter that is almost RS, except 
for a small perturbation of non-deltaic messages distributions: 

2 

Q[Q] = E p q s F (Q(u)-5 u , q ) (20) 

q=-2 

26 „ 

+ E / M{y})s F (Q{u)-Q n {{y})(u))d{y} 

n=l J 

and to see how it evolves. The symbol {y} stands for a subset of the parameters y2,yi,H-i, y~2 (see equation below). 
Following [28|, H^] this means that the order parameter has the following composition: 



2/2^,2 + (1 - V2)Su,l 
2/2<5„,2 + (1 - U2)Su,0 



with probability Pi,i e (— 2 . . . 2) 
with probability 711(2/2) 
with probability 7^(2/2) 



Q(u) = < y 2 5 Ut2 + yiS u ,i + (1 - Vi - 1/2)^,0 with probability wn(yi, Vi) 



(1 -t/i -y 2 - y-i - y-2)S u ,o+ with probability 7r 2 6(t/-2, V-i, Vi, yi) 

, y-i&u-i + y-2&u-2 



(21) 
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where the perturbation part is composed of all the 26 possible combinations of non-delta shaped distributions. We 
will refer to them as Q n {{y})(u) in general, and sometimes dropping the u for clarity. Q n {{y}){u) shall be interpreted 
as a function of u, where n and ({y}) are parameters, the former qualifying the type of non-deltaic distribution, and 
the latter defining the probabilities of each message in the distribution. Let us highlight that, while Q26({y}){ u ) 
seems to contain all the other possibles Q n ({y})(u), it is not the case, since, to be meaningful in the context of the 
decomposition (|21|) . all the components of the vector {y} must be greater than zero. This is a way to focus our 
attention directly on the overall weight of each type of perturbation: 

n„ = J ir n ({y})d{y} (22) 

For a first order study of the stability, the Tr n ({y}) are supposed to be infinitesimal quantities such that, during 
the iteration, the presence of more than one non-delta shaped distributions in the cavity will be very improbable. 
The convolution of a given Q n ({y}) with cK — 1 delta-shaped distributions in eq. (|12[) can yield, depending on 
the realization of the coupling constant, any of the distributions in (I2ip . Considering all the possibilities and their 
probability of occurrence in eq. (fT2"|) . we can construct the matrix 126x26 of elements i n>m giving the probability that 
a non-delta-shaped distribution of type n turns into type m after the iteration. Then, after the iteration, the vector 
IT = (III, Hj) • ■ ■ ! II26) of the overall weight of the perturbations is transformed by 

rr = in (23) 

A given RS solution is stable if all the eigen-values of I are smaller than one: 

max A(I) < 1 (24) 

This is the general approach in studying the RS stability. In some cases, this approach may lead to simple and closed 
equations that allow analytical solutions, but in general, the properties of the matrix I must be studied numerically. 



III. RESULTS FOR THE TRIANGULAR HUSIMI LATTICE 



In this section we apply the formalism described above to solve, at T — 0, the Edwards- Anderson model in the 
simplest of all Husimi lattices: the triangular lattice, with c = 2,K = 1 (see figQJ. Note that the case, p — —1 was 
recently extensively studied in [30j and is equivalent to the bi-coloring problem in an hyper-graph. Here we will keep 
K = 1 and concentrate our attention on the role of p in the thermodynamics of the system. 



A. RS-solution 



This graph allows closed expression for all the RS fixed points. For example, a careful analysis of (JT]) and (JSj) for 
this lattice leads to the definition of the following variables: mi = pi — p—\, Si = Pi + P-i, mi = pi — P-i and 
s 2 = P2 +P-2) such that Q may be written as: 

mi = (rnipo + y (mi + 2m 2 ))p - mi(y + s 2 ~Po)p 2 

, si. ml. m? 2 / s i\ s 

si = si[po+S2 + ^)-mi{m 2 -—)p+—p +{p -—)sip 

2 

P p 

m 2 = (siTOi + s 2 mi - sim 2 + 2m 2 )- + (miSi + s 2 mi + sim 2 + 2p rri2) — 

2p = s x + s 2 + s 1 s 2 + 2p + ni 1 m 2 p - (ml + m\ + m 1 m 2 )p 2 — (2po(l — Po) + «iS2)p 3 
s 2 = l-p - si (25) 

It is easy to check that, independently of p, a trivial paramagnetic solution (P) where only u — messages have a 
non-zero probability po = 1, always exists: 



Q P (u) = S U: q 



(26) 
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FIG. 2: The solution of the RS equations for the triangular Husimi lattice (c = 2 and K = 1). Panel A is the first spin-glass 
solution (SG-1). Panel B represents the ferromagnetic solution. In panel C appears the second spin-glass solution (SG-2). The 
variables are defined in section UlI Al 



On the other hand, setting the values of mi and si to zero, i.e. pi = p-\ = 0, automatically satisfies the first 
two equations in (|25|) . With the three remaining equations, we can find three solutions. The first one is a spin glass 
paramagnetic solution (SG-1) with po = (2p 3 + 1) , vti2 = 0, Si = 1 — pq: 

that is valid for p > psc~i — and is shown in panel A of figure [2J The second and third solutions are actually the 
same 2-degenerated ferromagnetic solution with po = s 2 = 1— Po and m| = (p 2 +p— l) 2 p~ 6 — 2(p 2 +p— 1)(1— p)p~ 3 '- 

s2 - to2 1 - p c s2 + m2 

Qf " = ^ o u ,-a H 5— <Vo H ^ U 2 (28) 

2 /cr 2 

valid for p> pp = ~ 0.7071. The dependencies of po 5 s 2 and m_2 with p appear in panel B of figure^ 

There is yet another spin glass paramagnetic solution with s\ ^ 0, toi = and vri2 = that we will call SG-2. It 
differs from the SG-1 in that p\ = p-\ > 0. The dependencies of si, S2 and po for this solution are shown in the panel 
C of figure [51 Like the other solutions described above, this solution can be analytical expressed as a function of p. 
However, the resulting algebraic expression is quite large and does not add clarity to the remaining discussion, so we 
do not write it explicitly here. Just note (figure [2]) that this solution exists only for p > psc~2 = 0. At this point it 
is worth mentioning that the notations SG-1 and SG-2 refer to the non-trivial structure (although symmetric) of the 
distributions, but one must keep in mind that the solution is ergodic and does not represent a true glassy system. 

Each of the four solutions described above defines an energy function U (p) by eq. (see figure \5§ . Except for the 
second spin glass paramagnetic solution, which we leave to the interested reader to work out, the expressions for the 
energy are rather simple: 

Up = -^(2 + p 3 ) (29) 
2 2 + 9p 3 + 12p 6 

USG-1 — -77- 



3 (i + 2 P 3; 



2 



2 (-1 + 3p - 7p 3 + 6/9 4 - p 6 + 3p 7 ) 

U F = —-- 



3 3p 6 

Summarizing, the RS calculations show that: i) a trivial paramagnetic solution exists for all p; ii) at p — psc~i = 
Psg-2 = two spin glass paramagnetic solutions appear and are valid up to p = 1; Hi) for p > pp = 2 -1 / 2 one also 
finds a ferromagnetic solution. Note that, in the scale of the figure [31 it is hard to differentiate the actual energies of 
the different solutions for a wide range of values of p. As we demonstrate in the next section, the same is true when 
compared with the 1RSB solution. However, we will show that for other lattices the situation is different. 
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FIG. 3: Energies for the different RS solutions and the 1RSB approximation 4>(fJ,*) as a function of p for the triangular Husimi 
lattice (c = 2, K = 1). We represent also the stability threshold pp = of the trivial paramagnetic solution. The caption 
shows the point in which the ferromagnetic solution appears pF, and the point where it becomes stable p%. 



B. Stability analysis of the RS solutions 



Though the general method given in the introduction is always valid, we will try to keep the stability analysis as 
simple as possible. The trivial paramagnetic solution can be studied, for instance, with only one general perturbation 

Q26(y-2,y-i, 2/1,2/2) = V-2&U-2 + y-iSu-i + (1 - y-2 - y-i -yi- 2/2)^,0 + yiS u ,i + 2/2^,2 

Q[Q] = S F {Q(u)~S(u,0)) (30) 
+ J K(y 2 ,yi,y-i,y-2)5 F (Q(u) - Q26(y2,yi,y-i,y-2))dy 2 dyi . . . 

The idea is to follow the overall weight of the perturbation II = J 7r(j/2, J/ij V-i, 2/-2)dj/2dj/idj/_idj/_2 for the iteration 
(fTSjl . The reason behind the success of such a simple analysis lies in the fact that the iteration eq. p^|) has a very 
naive behavior when we consider the convolution of the paramagnetic solution with a perturbation. Depending on 
the realization of the coupling constants Jij, the convolution of the perturbation Q2&{y-2, y—i, yi, 2/2) and the RS 
solution S(u, 0) can only produce one of the following three results: 

5 U ,0 * Q26(2/-2,2/-l,2/l,2/2) -> S ufi 

5 u ,o*Q26(y-2,y-i,yi,y2) -» 026(2/2, yi, y~ uy-2) 
6 Ut Q*Q 2 6(y-2,y-i,yi,y2) -> Q 2 6(y-2,y-i,yi,y 2 ) 

where we are using the symbol * to represent the right-hand side of eq. (|12p . Averaging over all possible couplings, we 

found that the last two situations take place with a probability . Then, also taking into account a combinatorial 
factor of 2, the overall weight of the perturbation evolves following: 



u'p = (1 + p 3 )ii P (31) 

showing that for all — 1 < p < pp = the paramagnetic solution Qp{u) — 8 U: q is stable, while it is unstable for p > 0. 
To determine the stability of the solution at pp — one must use a second order perturbation. 
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Unfortunately, in general, such a simple analysis is not applicable. For instance, the spin glass paramagnetic 
solution, SG-1, when convolved with Q26! can yield mixed distributions of other types. Then, in order to keep 
analytical expressions, we consider in this case two other types of perturbations: 

Qn(y) = (l - y)S u ,a + y$u,2 

Qm{y) = ySu-2 + (1 - y)5 u ,o (32) 

Both the spin glass paramagnetic solution, SG-1, and the ferromagnetic solution, F, have only messages u 6 
{—2,0,2}. Then, the convolution with the perturbations Q n ,m produces one of the following results: 



QSG- 


-1,F 


* 




-* 8 u fi 


QsG- 




* 


Qn,m(y) ~ 


S u< 2 


QsG- 




* 


Qn,m(y) ~ 




QsG- 


-1,F 


* 


Qn,m(y) ~ 


■+ Qm(y') 


QsG- 


-1,F 


* 




+ Qn(y') 



where only the last two situations are responsible for the evolution of the perturbation. This set of perturbations is 
closed under the iteration equation 1|12[) and therefore suitable for analytical treatment. 

For the spin glass solution, SG-1, the matrix for the evolution of the perturbation lf^ 2 ~ 1 is symmetric, with elements: 



_ 1 + 3p 3 + p + p 2 - 2p 5 

2(1 + 2p 2 ) 
1 + 3/r 3 - {p + p 2 - 2p 5 ) 



2(1 + 2p 2 ) 

l+2p :i 



The largest of its eigenvalues is A = 1+9^ which is always greater than one, for p > 0. This means that the spin 



glass SG-1 solution is always unstable. 

For the ferromagnetic solution, the matrix lf x2 has the following structure 



3p^- 1 I 
2p 



1 -(l + --p-2p 2 ) 



itm* Vn = ± -^V(P 4 -(P-1) 2 )(V-1) 



where the terms *^" m ,* m „ differ in the sign in front of the square root. The limiting condition for the stability 
max A (I) = 1 turns out to be equivalent to the solution of p 6 + p A — 2p 3 + p 2 + p — 1 = 0. This defines the point 
Pp ~ 0.765942 where the ferromagnetic solution becomes stable. 

For the spin glass paramagnetic solution SG-2 it was impossible to find a closed subset of perturbations smaller 
than the full set. Therefore we studied its stability using the full matrix 126x26 ■ Given the size of the matrix, and 
the non-simple dependence of this solution with p, the study of the stability was done numerically. We calculated all 
the 26 eigenvalues of 126x26 and found that some of them were always above 1 in the interval < p < 1. This means 
that the spin glass solution (SG-2), like the (SG-1), is always unstable. 

Note, that the results obtained for the SG-1 and ferromagnetic solutions where only a subset of perturbations were 
studied, can be considered only as lower-bound approximations for the stability. For the SG-1 the full 26 perturbations 
procedure can be done analytically, and it appears that the highest eigenvalue was, in fact, the one we found with 
the restricted perturbation method. For the ferromagnetic solution it is harder to get an analytic expression for the 
eigenvalues of the full matrix 126x26: but yet we checked numerically that all the eigenvalues of this matrix were indeed 
equal or smaller than the one found within the restricted perturbation analysis. In summary, the triangular Husimi 
lattice with K = 1 has a stable RS trivial paramagnetic solution for —1 < p < pp = 0, has a stable ferromagnetic 
solution Qp(u) when pp < p < 1, and two spin glass paramagnetic solutions that are unstable to replica symmetry 
breaking. These points of stability are shown in figure 



C. 1RSB solution 



The 1RSB approximation for the free energy of a model is a formalism more general than the Replica Symmetric 
calculation. Thus, in the calculation of the free energy, we do not need to restrict ourselves to the interval pp < p < pp 



11 



where the RS solutions are unstable. Instead we can do the 1RSB calculation in the whole interval — 1 < p < 1, 
and the RS solutions will still be found if they are thermodynamically significant, which is a stronger condition than 
stability. 

Using the population dynamics described in section fll A 2l we computed 4>{p) for different p. In figure[3]we compare 
the ground-state predictions Urs and Uirsb for both approximations (RS and 1RSB) as a function of p. In the 
caption we show how the 1RSB solution collapses into the F solution for p > pp. Unfortunately, with the resolution 
of figure [3] it is impossible to differentiate the ground state energy in the 1RSB approximation from the energy of the 
RS solutions within the region in which no RS solution is stable, although the order parameter Q[Q] showed to be 
non-trivial, i.e. not having the form of (119)) . 

Triangle c=2 K=l 

-1.240 

-1.250 

-1.260 

-1.270 
-1.49 
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"°" -1.492 

-1.493 
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-1.820 
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H 



FIG. 4: Free energy <j)(p) of the 1RSB solution for the triangular Husimi lattice (c = 2, K = 1) for different values of p. 
Top panel: paramagnetic region, Middle: 1RSB region and Bottom: ferromagnetic. The flat curves are representative of RS 
solutions, while the curve in the 1RSB region has a clear non trivial maximum 4>(p*) = Uirsb at /i ~ 0.6. 

For the sake of understanding the qualitative differences between the RS and the 1RSB solution, we show the curve 
</>(/i) for three different p in figure 2] The lower curve corresponds to the ferromagnetic region (p > -^), the upper 

to the paramagnetic (P) interval, and the middle curve to the 1RSB zone. In the paramagnetic and ferromagnetic 
regions, 4>(p) is flat, as expected. It has zero derivative and hence zero complexity, meaning that the solution space is 
not clustered and that the RS assumptions are valid. Furthermore, a detailed analysis of the order parameter shows 
that it actually has the RS structure (fl9|) . On the other hand, in the 1RSB region (p = 0.65) 4>{p) is convex and has a 
maximum at p = p* ~ 0.6, and the order parameter is genuine 1RSB (not made of deltaic distributions) for all values 
of p. The highest point of this curve defines the 1RSB approximation for the ground state energy Uirsb = 0(/ i *)- 

The learned reader will find rare that a highly-frustrated system, as the Husimi triangular lattice with an excess of 
antiferromagnetic interactions (p < 0), does not have a 1RSB solution. He might argue that the proved stability of 
the trivial paramagnetic RS solution is not enough to rule out other solutions. In fact, the absence of 1RSB solution 
may well be an artifact due to the random initial conditions used in the population dynamics. Fortunately, with the 
same tools used to analyze the stability of RS solutions, it is possible to prove (see appendix A) that any Q[Q(u)] 
solution of the self-consistent equation (jlip is unstable to a paramagnetic perturbation Q[Q] = ir8 F (Q{u) — S Ui o). 
This rules out any possible stable solution except for the complete paramagnetic one Qp[Q(u)] — 5 F (Q(u) — <5 u ,o), 
rendering the trivial paramagnetic solution the only thermodynamically relevant in the region — 1 < p < 0. It must 
be kept in mind, however, that at finite temperature a 1RSB solution with real-valued fields, which are not embedded 
in our parameterization, may facilitate a RS-1RSB transition in this zone. 

To help the intuition, let us look closer to the case p = — 1. In this case, each triangle has a minimum energy 
6-times degenerated, and only 2 excited states. It is true that in all the states of minimum energies, one of the bonds 
of the triangle is violated, but one cannot forget that this is anyway the ground state of the triangle. Then, because 
this state is highly-degenerated, it is not hard to imagine that there is a lot of freedom, even on the large scale of the 
network, to have all the triangles "satisfied" . Of course, this argument is valid provided K is low. If K is high enough, 
each spin will be connected to a large number of triangles, and therefore, the freedom of choice will be reduced and 
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eventually may disappear. This, in fact, should turn the problem to the non-trivial 1RSB scenario predicted in[30| 
for K > 6. 

IV. GENERAL HUSIMI LATTICES 

The procedure described in the introduction, and applied to the study of the triangular Husimi graph, is also 
applicable to all other Husimi graphs. The main difference arises in the difficulty in generalizing the analytical results 
obtained above. For instance, the equation (|7|) is still parameterized by P-27P-I7P07P1JP21 but since each term has 
degree c • K on these probabilities, it is very hard to write closed equations in terms of p. The absence of closed 
expressions for these probabilities also makes it impossible to express the RS-1RSB instability analytically. However 
we can always solve numerically the set of five equations (O on the probabilities pi for each value of p, and check 
numerically the eigenvalues of the stability matrix for such solutions. The only exception that allows a full analytical 
treatment independently of K and c, is the trivial paramagnetic solution Qp = S U fl, that is obviously a solution of 
equation for all the Husimi lattices. 

A. The trivial paramagnetic solution 

The result obtained for the trivial paramagnetic solution in the cactus can be easily generalized to all kinds of pure 
and regular Husimi lattices, with generic loop size c + 1 and degree K + 1. The paramagnetic solution is always valid, 
but it changes from one type of graph to another. To study the transition RS-1RSB we consider again a perturbed 
order parameter of type (|3T))) . In the general case equation (fT2")l is a convolution of c • K distributions. Again, the 
convolution of paramagnetic distributions reproduces a paramagnetic distribution: 

($0,u*) CK -> V« 

The non-trivial case is the convolution with the perturbed distributions Q26 that are present in the order parameter 
Q[Q] in an amount proportional to n2e{y2, ViiV-i, y-2) (see equation (|30|) ) . Depending on the realization of the 
coupling constants, the integration of eq. Q12[) may have three different results: 

(<5«,o*) c ~ x Q26(y2,yi,y-i,y-2) -> S Ul0 

(Sufi*)^^ 1 Q2e(y2,yi,y-i,y-2) -> Q 2 6(j/2,yi,y-i,y-2) 

(5u,o*) c Q26(y2,yi,y-i,y-2) -> Q2e(j/-a,j/-i, 2/1,2/2) (33) 

It is easy to check that, independently of the loop size c + 1 and the connectivity K + 1 of the Husimi graph, 
the last two outcomes (those that are responsible for the propagation of the perturbation) take place only when an 
even number of couplings Jij are anti-ferromagnetic. The probability to have an even number of anti-ferromagnetic 
couplings around a loop can be computed as: 



(c+l)/2 . . c+l-2q /, \2q 



1 + pV l -p\ ( 2 i \ _ 1 



o(l + P C+1 ) (34) 



^ \ 2 \ 2 \ c+1 2 

q=0 

The upper limit of the sum is taken to be the higher integer value not greater than (c+ 1) /2. Since the perturbation 
Q26 can occur in any of the c- K cavity biases distribution of eq. (|12p .one needs to consider a cK combinatorial factor 
in the previous expression. Thus we get that the perturbation weight Hp = J 7r(j/2, J/i , y~i, J/-2) evolves through: 

rK 

n P = — (i + P c+1 )n P (35) 

The condition of stability is the damping of the perturbation ^j-(l + p c+1 ) < 1 which has no solution for odd c, i.e. 
for even loop sizes as squares, hexagons, etc. The general solution: 

s (cK-2\^ 



P<Pr = -[^ !r ) (36) 

is valid for odd loop sizes like the triangle, the pentagon, etc. For values of p < pf, the paramagnetic solution 
Q(u) — Sq.u is stable. On the other hand, for even loop sizes the paramagnetic solution is never stable. 
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FIG. 5: Energies for the different RS solutions and the 1RSB solution as a function of p for the square Husimi lattice (c = 3, 
K = 1). The appearance and the stability of the solution F are signaled. The caption shows the difference between the 1RSB 
approximation Uirsb ~ 0(/Lt*) (black squares) and the energy of the RS solutions. White squares represent 0(0) and coincide 
with the energy of the SG-2 solution 



The energy of the paramagnetic solution can be computed very easily. Given that only u = messages exist, each 
cavity loop preserves the symmetry of the Hamiltonian (S <-> —5), and then the only contribution to the energy 
comes from the frustration inside the loops. Frustrated loops are those with an odd number of anti-ferromagnetic 
interactions, and have an energy — (c — 1), while loops with an even number of anti- ferromagnetic interactions are not 
frustrated and have energy — (c + 1). Then, using (|34[) we get as the energy of the paramagnetic solution: 

U P = -^±±(c + p c+1 ) (37) 
c + 1 

So, two key difference arise between Husimi lattices with even and odd loop sizes. The parity of the energy 
function U p (p) is exactly that of the loop sizes c + 1 , i.e. graphs with even loop sizes have an even energy function 
Up(p) = Up(—p), while graphs with odd loop sizes have an odd energy dependence with p (except for a constant), as 
shown in figures [5] and [6] Furthermore, the trivial paramagnetic solution is stable (under 1RSB perturbations) only 
in lattices with odd loop sizes. 



B. Spin glass and ferromagnetic solutions 

The numerical solutions of {7J for different values of c and K show a zoology of phases similar to the one found 
in the triangular lattice: a trivial paramagnetic solution (P) already discussed in detail, two spin glass paramagnetic 
solutions (SG-1) and (SG-2), and a ferromagnetic solution (F) (see figures [S] and O . The solutions SG-1 and F give 
zero probability to messages of type u = ±1, as in the triangular Husimi lattice. 

As expected, for a sufficiently ferromagnetic systems, p > pf, there is a ferromagnetic phase in all types of Husimi 
graphs. The spin glass solutions, however, are related to the parity of the loops. Husimi lattices with even loop sizes 
(square c+1 — 4, hexagon c+1 = 6, etc ) present the spin glass solutions for all p (see figure [SJ. On the contrary, in 
odd loop sized Husimi graphs, the appearance of the spin glass solutions seems to be related to the loss of stability of 
the trivial paramagnetic solution (see figure [6]). Other Husimi lattices with odd and even number of loop sizes c+1 
and K = 1, 2 reproduced the same behavior. 

The stability of these solutions was analyzed numerically by studying the eigenvalues of the 126x26 matrix (see section 
III Bp . In general the panorama is similar to the one found for the triangular case; the two spin glass paramagnetic 
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FIG. 6: Energies for the different RS solutions and the 1RSB solution as a function of p for pentagonal Husimi lattice (c = 4, 
K = 1). The stability threshold for the P and F solutions are signaled. The caption shows difference between the 1RSB 
approximation Uirsb = 4>(p*) (black squares) and the energy of the RS solutions. White squares represent (f>(0) and coincide 
with the energy of the SG-2 solution. 

solutions are always unstable to 1RSB perturbations, and the ferromagnetic solution has a non-trivial stability point 
pp above which it becomes stable (sec figures [5] and 15]) . In table (JTJ) we present, for different c and K, the numerical 
values for these stability points and for the appearance of the ferromagnetic solution pp. 



Husimi 


Pp 




4 


c= 2, K = 1 





1/V2 ~ 0.707 


0.766 


c = 3, K = 1 




0.671 


0.749 


c = 4, K = 1 




\[\ ~ -0.8706 


0.665 


0.748 


c = 2, K = 2 


— -1 


^1 ~ -0.7937 


0.534 


0.631 



TABLE I: The paramagnetic and ferromagnetic solutions are stable to the left and right of pf> and pp respectively. The 
stability of the paramagnetic solution is found analytically (see section Hi Bp . The point in which the ferromagnetic solutions 
first appears within the RS approximation pp, as well as the stability of this solution, are found numerically by solving the RS 
self consistent equation |7]) for different p. 

The picture that emerges from these calculations is the following: above pp, the ferromagnetic RS solution is stable 
to 1RSB perturbations and thermodynamically relevant. Below this point, both spin-glass solutions are unstable. On 
the other hand, the stability of the trivial paramagnetic solutions (P) depends on the parity of the loops. In Husimi 
lattices with even loop sizes (c+ 1 = 4, c+ 1 — 6, etc) this phase is never stable, while in lattices with odd loops sizes, 
this phase becomes stable below pp. 

C. 1RSB Solutions 

For each Husimi lattice and for each value of p, the 1RSB free energy <p(pi) is computed using the usual method 
of population dynamics. The 1RSB prediction for the energy of the model is given by Uirsb — 4>{^*) where pf 
maximizes <f>. On the other hand, the limit p, — ► must recover the RS equations and 4>(0) must coincide with one of 
the replica symmetric solutions. 
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In figures [5] and [6j Uirsb — </>(/•**) and 0(0) are shown simultaneously for the square and pentagon Husimi lattices 
with K = 1. In both types of graphs, the ferromagnetic solution is thermodynamically relevant for large p, as our 
intuition tells us. In the case of odd loop sizes, the situation is similar to the one found in the triangular Husimi 
lattice: the trivial paramagnetic solution seems to be the thermodynamically relevant for negative enough p. 

IRSB-RS-Ferro Transition in c=4 K=l 
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P 



FIG. 7: Energy versus p close to the IRSB-Ferromagnetic transition. The real thermodynamic behavior of the system is defined 
by the black squares. The replica symmetric ferromagnetic phase is not only stable above p%, but also thermodynamically 
relevant, as the 1RSB solution smoothly becomes RS (see eq. I19[) as p — > p%. 

Summarizing, below the ferromagnetic stable region pf, < p < 1, the model has always a thermodynamically 
relevant 1RSB solution. In Husimi lattices with even loop sizes (c+ 1 = 4, c+ 1 = 6, etc.) this phase goes till p = — 1, 
while in Husimi graphs with odd loop sizes and depending on K, this 1RSB phase may disappear for negative enough 
p and in this case, the thermodynamics of the system is controlled by a trivial paramagnetic solution. However, it is 
worth mentioning that for large cK, pp — » —1. Moreover, we will show below that the 1RSB solution may extend 
below pp, so, in more general cases, we expect that this solution becomes thermodynamically relevant also for p = — 1 

& 

To understand the characteristics of the transitions between the RS to the 1RSB phases we take the c = 4, K = 1 
case as a model. We made a close-up and a careful calculation around the IRSB-Ferromagnetic and the 1RSB- 
Paramagnetic transitions. They appear in figures [7J and [5J respectively. In both cases the simulations were done 
adiabatically, starting from the 1RSB zone and slowly varying p towards the transition points and working always at 
p* . In this way we guarantee that the algorithm finds (if it exits) the 1RSB solution. 

As the figure [7] suggests, the IRSB-Ferromagnetic transition is continuous. Increasing p, a close inspection of the 
structure of the order parameter Q[Q(u)] (represented by the population of distributions Q(u)) shows that the amount 
of deltaic distributions Q(u) = 5 u ,u grows smoothly as p approaches pp. 

On the other hand (see figure [8j, the 1RSB- Paramagnetic transition is quite different. There is region pmsB < 
p < pf, in which the thermodynamically relevant solution is 1RSB, but where the RS solution is stable under 1RSB 
perturbations. Then, below the point Pirsb ~ —0.89 the paramagnetic solution has higher energy and becomes 
thermodynamically relevant. Note on the inset that below this point there is still a region where the algorithm 
finds a non-trivial 1RSB solution of lower energy. This suggests that in this zone 1RSB solutions may appear but 
are exponentially rare. Then, by further decreasing p, the population dynamics finds only the trivial paramagnetic 
solution. The inspection of the distribution shows that this transition is discontinuous. Above Pirsb the order 
parameter Qirsb[Q(u)] is free of deltaic distributions Q{u) — S UtUo and below Pirsb the solution is completely 
paramagnetic. 

D. The limit of the Bethe lattice 

Finally, we focused on the unbiased model p — 0. In this case, all the Husimi lattices, with the exception of the 
triangular (c = 2 and K = 1), are at least 1RSB (see eq. V 
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RS-1RSB Transition in c=4 K=l 
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FIG. 8: Energy versus p close to the IRSB-Paramagnetic transition. The real thermodynamic behavior of the system is 
defined by the black squares, except for the unstable region in which the adiabatic continuation of the 1RSB solution goes 
below the paramagnetic curve (see caption). For p < pirsb the thermodynamically relevant order parameter is the RS trivial 
paramagnetic solution. The transition IRSB-Paramagnetic is not continuous: the structure of the order parameter changes 
abruptly when crossing pirsb- 



If wc refer to the Husimi lattice considering only the local connectivity of the sites (number of neighbors), then, a 
Husimi lattice with degree K ± 1 is locally very similar to a Bethe lattice with degree k + 1 = 2(K + 1). In table ITT1 
we present a few ground state energies at the 1RSB level of approximation for different Husimi lattices and compare 
them with similar results for Bethe lattices. Note that for c = 7, K = 1 and c = 5, K = 2, the differences between the 
Bethe and the Husimi ground states are already very small. 



Husimi 


Ugs 


c=2, K =1 


-| = -1.333 


c = 3, K=l 


-1.444 ±0.002 


c=5, K=l 


-1.468 ±0.002 


c = l,K=l 


-1.470 ±0.002 


Bethe, k — 3, c = oo,K = 1 


-1.471 ±0.002 



Husimi 


Ugs 


c = 2, K = 2 


-1.471 ±0.002 


c = 3, K = 2 


-1.819 ±0.002 


c = 5, K = 2 


-1.823 ±0.002 


Bethe, k = 5, c = oo,K — 2 


-1.825 ±0.002 



TABLE II: Ground state energies within the 1RSB approximation of different Husimi lattices for p — 0. 



This similarity goes beyond the ground state calculations. In figure [9] we plot the free energy </>(//), and the 
complexity E(e) for different Husimi lattices with K = 2 and p — and for the Bethe lattice. As expected, the figure 
shows a clear maximum at <f>(p,*) and the two usual branches in the S(e) plot[12j. From these plots, we may also 
conclude that the local similarity between Bethe and Husimi lattices, for which the relation k + 1 = 2(K ± 1) holds, 
turns out to be true, on thermodynamic grounds for loop sizes greater than c±l = 6(c±l = 8if-ftT = l). This 
suggest that, at least for the EA model, the Bethe approximation may work even under conditions less restrictive 
than loop sizes of order \nN. 



V. CONCLUSIONS 



We solved the Edwards- Anderson model in a Husimi graph at T — 0. We presented closed analytical (RS) expres- 
sions as a function of p for the existence and the stability of a trivial paramagnetic solution (P). For the triangular 
Husimi lattice (c = 2 and K — 1) we obtained similar expressions for the appearance and the stability of two spin-glass 
solutions (SG-1, SG-2) and a ferromagnetic solution (F). For other cases, these points were calculated numerically 
(see table U). Within the 1RSB approximation we obtained, using a population dynamics algorithm, the value of the 
ground state energies of different lattices as a function of p (see figures [5] and . 
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FIG. 9: Free energy (f)(p) (left) and complexity E(e) (right) of different Husimi lattices at p = and K — 2. The graphics for 
K — 1 axe similar. Also represented the EA solution in the Bethe lattice with degree k — 5. The symbols in the right panel 
correspond to the legend in the left panel. The early coincidence (c = 5) between our short loop model in Husimi graphs, and 
the short-loop-free model in the Bethe lattices, gives a thermodynamic threshold for the relevance of loops. 



The main picture emerging from this work is that for Husimi lattices with even loop sizes, the system is at least 
1RSB in all the range — 1 < p < pp(c,K). Above p F (c,K) the model has a RS ferromagnetic solution, and the 

transition from the 1RSB phase to the ferromagnetic one is continuous. On the other hand, for lattices with odd loop 

i 

sizes a trivial paramagnetic solution is stable under 1RSB perturbations up to pf, = — ( y 2 ) c+1 , in addition for 
small K this solution may become thermodynamically relevant below Pirsb < Pp where the order parameter jumps 
discontinuously from a non-trivial 1RSB solution towards trivial paramagnetic distributions. For the particular case 
of the triangular lattice with K = 1 it was proven that the trivial paramagnetic solution is actually the unique and 
thermodynamically relevant solution (at least in the 1RSB frame) for all p < 0. 

Finally, we focused on the p = case and computed the ground state energies for different values of c and K (see 
table |TT| . Our results suggest that the energy and complexity of the EA model in Husimi lattices with loop sizes 
c + 1 > 8 are already well described by an equivalent short-loop-free EA model in a Bethe lattice. This, in turn, can 
be considered a thermodynamic threshold for the shortness of loops, less restrictive than the usual ~ \og{N). 
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VII. APPENDIX A 



Let us prove that the trivial paramagnetic solution is actually the only thermodynamically relevant solution of the 
self-consistent equation (|11[) for the triangular Husimi lattice (c = 2 and K = 1), in the interval (—1 < p < pf, = 0). 
In doing so, we will follow exactly the same procedure used to study the instability of the Replica Symmetric solutions. 
We will presume that a non-trivial solution exists for the self-consistent equation, and we will show that it is unstable 
to small variations of the probability of paramagnetic distributions of messages Q{u) = <5 u .o. 

Let us call Q[Q] the solution of the self-consistent equation at 1RSB level of approximation. A perturbed order 
parameter would be Q[Q] + ttS f (Q(u) — S u .o). If we write down the self-consistent equation (fTTj) using this order 
parameter, and keep to the first order in it in the right-hand side, we will get the following equation for the evolution 



18 



of the perturbation weight: 

7r'<Su,o = 27rE,/ ,o,Q,du,o])dQ[Q] (38) 

where Ej is the expectation over the coupling constants J, and the 2 multiplying the integral is a combinatorial factor. 
The possible realizations of the disorder can be grouped into two sets: those with an odd number of antiferromagnetic 
interactions, like (—1, —1, —1) and (—1, 1, 1); and those with an even number like (1, 1, 1) and (—1, —1, 1). Elements 
of the first group occur with a probability 

^frustration - l~ ^ > + 6( -^~ ' ' ~ "IT" ^ ' 

while the other group has a complementary probability. Whenever the disorder happens to be frustrated (first group) , 
the convolution Q[J_, Q, S Ui o] = <5„,o regardless of the actual distribution Q. On the other hand, if the disorder around 
the triangle is not frustrated (second group) Q[J_, Q, <5 Uj o] = Q(u), thus reproducing the probability distribution of 
Q{Q{u)\. Then the contribution to the perturbation is only given by the frustrated triangles, and we get for the 
iteration: 

tt' = 27rP frustraci6n (40) 

The instability condition would be given by: 

^frustration — ^ — P > ^ 

which happens to be true in the interval — 1 < p < 0. This proves that any order parameter is instable to paramagnetic 
perturbations. This leaves no room for other stable order parameters than the trivial paramagnetic one, concluding 
our proof. 
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